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The bare and hydrogen-covered diamond (100) surfaces were investigated through pseudopo- 
tential density-functional calculations within the local-density approximation. Different coverages, 
ranging from one to two, were considered. These corresponded to different structures including 1x1, 
2x1, and 3x1, and different hydrogen-carbon arrangements including monohydride, dihydride, and 
configurations in between. Assuming the system was in equilibrium with a hydrogen reservoir, the 
formation energy of each phase was expressed as a function of hydrogen chemical potential. As the 
chemical potential increased, the stable phase successively changed from bare 2x1 to (2xl):H, to 
(3xl):1.33H, and finally to the canted (lxl):2H. Setting the chemical potential at the energy per 
hydrogen in H2 and in a free atom gave the (3xl):1.33H and the canted (lxl):2H phase as the 
most stable one, respectively. However, after comparing with the formation energy of CH4, only the 
(2xl):H and (3xl):1.33H phases were stable against spontaneous formation of CH4. The former 
existed over a chemical potential range ten times larger than the latter, which may explain why the 
latter, despite of having a low energy, has not been observed so far. Finally, the vibrational energies 
of the C— H stretch mode were calculated for the (2xl):H phase. 

PACS numbers: 71.15.Nc, 82. 65. My, 82.20.Wt, 68.35.Md 
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I. INTRODUCTION 



Recently the nature of the structures of hydrogen on 
diamond (100) surfaces has attracted new interest due to 
technological advances of epitaxial growth of this surface 
by chemical vapor deposition (CVD) methods. In the 
CVD homoepitaxial growth of diamond, the surface is 
mostly hydrogen terminated because abundant hydrogen 
exists in the growth environment.^ Atomically smooth 
surfaces were obtained only in the (100) orientation.'*^® 
Therefore, understanding the structural and electronic 
properties of this surface becomes quite important. Since 
there are two dangling bonds per surface carbon atom for 
the bulk-truncated bare surface, a reconstruction takes 
place in which two adjacent surface atoms move closer 
toward each other to form rows of symmetric dimers. 
Upon adsorption of hydrogen atoms, possible configura- 
tions include the monohydride (one hydrogen attached to 
a carbon atom) and dihydride (two hydrogen atoms at- 
tached to the same carbon atom) arrangements. The for- 
mer is expected to keep the carbon dimer reconstruction 
(2x1), while the latter saturates both dangling bonds of 
a surface carbon atom, yielding a structure likely to be 
1x1 and high in energy due to the steric repulsion. 

Unlike the bare surface, what really happens on the 
hydrogen-covered C(IOO) surfaces was under debate for 
many years. A definite experimental determination of the 
stable structure was hindered by the difficulty of produc- 
ing smooth surfaces, i.e. surfaces with large, atomically 
flat terraces. ^^""^"^ Only recently have researchers resolved 
this problem.'*^® Thomas and Butler** reported high res- 
olution electron energy loss (HREEL) spectra indicating 
only monohydride species existed after ex situ exposure 



of the surface at 800° C to a hydrogen plasma and no di- 
hydride was observed after m situ exposure of the surface 
to hydrogen atoms produced by a hot tungsten filament. 
In addition, they observed a sharp two-domain 2x1 re- 
construction by low energy electron diffraction (LEED) 
for the hydrogen saturated surface. Kiittel et al.^' ana- 
lyzed natural, doped (100) diamond surfaces polished in 
a hydrogen plasma by atomic force microscopy (AFM), 
LEED, and X-ray photoelectron diffraction. They found 
that after etching the C(IOO) surfaces had a roughness 
of 0.8 nm and showed a sharp 2x1 LEED pattern which 
was even stable in air without annealing, and that ex- 
posing these surfaces to atomic hydrogen produced by 
a heated filament led to an increase of surface roughness 
comparable to what was observed on the as-received sam- 
ple. They concluded that the often reported 1x1 LEED 
pattern seemed to be a consequence of a large surface 
roughness and was absent on smooth surfaces. 

Using the vibration sum-frequency generation (VSFG) 
spectroscopy, Anzai et al}'^ observed vibrational reso- 
nance peaks at 2910 and 2960 cm~*. They assigned the 
main 2910 cm~* band to the CH stretch vibration of 
the monohydride groups on the C(100)(2xl) surface, and 
the 2960 cm~* band to the anti-symmetric CH stretch- 
ing vibration of the dihydride groups located at either 
the unreconstructed C(IOO) (1x1) island surface or step 
sites. Kawarada et al}^ examined the 2x1/1x2 surface 
reconstruction of C(OOl) surface using a scanning tunnel- 
ing microscope (STM) at an atomic scale and reflection 
electron microscopy (REM) at a microscopic scale. 

Quite a few theoretical studies of the structural and 
electronic properties of the H-covered C(IOO) surface 
have been performed, ^^^^^ including two recently pub- 
lished first-principles slab calculations**'*® based on the 
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local-density approximation (LDA). Furthmiiller et al}"^ 
used the ultrasoft pseudopotentials^'* and the conjugate- 
gradient technique, while Zhang et al}'' employed the 
Car-Parrinello molecular-dynamics (MD) method. Both 
examined the (2xl):H monohydride structure and found 
it to have a large desorption energy. Furthmiiller et 
al}"^ also calculated a higher-coverage structure with 1.5 
monolayer (ML, with respect to surface carbon atoms) of 
hydrogen on the surface, and found that the desorption 
energy for it to go to the (2xl):H monohydride phase 
was rather small. Zhang et al}'' reported that the dihy- 
dride phase is energetically favored over the monohydride 
phase relative to H desorption, while it is unfavorable 
with respect to H2 desorption. On the other hand, ear- 
lier studies by Yang et al}^' using a nonself-consistent 
tight binding model with parameters obtained from den- 
sity functional theory (DFT), and by Skokov et al}'' us- 
ing molecular dynamics with a semiempirical Hamilto- 
nian, reported a (3xl):1.33H structure to be most favor- 
able energetically. This structure consists of alternating 
monohydride carbon dimers and dihydride units. It was 
previously observed^®'^® in the hydrogen-covered Si(lOO) 
surface, but has not been identified experimentally for 
C(IOO). A fully self-consistent calculation is therefore 
necessary to clarify this situation. In addition, in or- 
der to compare stability of different phases with differ- 
ent hydrogen coverage, one has to take into account the 
chemical potential of hydrogen. In this study, we will 
evaluate the energies of different structures with cover- 
ages [9) ranging from one to two ML, including all the 
phases mentioned above. Depending on the values of the 
chemical potential, the equilibrium phase in contact with 
a hydrogen reservoir can be the (2xl):II, (3x 1):1.33II, 
or (lxl):2II (canted) phase. Comparing these forma- 
tion energies with that of a stable hydrocarbon CII4, one 
can exclude the (lxl):2II phase. Then the stable surface 
phase is (3x 1):1.33II for a small chemical potential range 
and (2x 1):II for a wide chemical potential range. We will 
also present the energetics in terms of the desorption en- 
ergy and compare the results with previous studies. Fi- 
nally, the vibrational energy of the C— H stretch mode is 
calculated for the (2xl):II phase. 



II. CALCULATIONS 

The calculations were carried out within the local- 
density approximation (LDA)^® with the Ceperley-Adler 
exchange-correlation potentiaP^ parametrized by Perdew 
and Zunger.^*^ The soft pseudopotentials were generated 
by the scheme of Troullier and Martins^^ for both car- 
bon and hydrogen. The carbon pseudopotentials were 
generated in the ground state configuration 2s^2p^ . The 
cutoff radii were 1.50 and 1.54 a.u., respectively, for s and 
p potentials. It has been shown that good convergence 
can be achieved with a plane-wave cutoff of 49 Ry for 
these pseudopotentials.^^ The hydrogen pseudopotential 



was generated from the Is^ configuration with a cutoff 
radius of 1.25 a.u. In the Troullier and Martin scheme, 
the first four derivatives of the pseudo- wave-function are 
matched to those of the all-electron wave function at the 
cutoff radius. Hence noticeable deviations of the pseudo- 
wave-function from the all-electron wave function nor- 
mally occur at positions a few tenths of an a.u. inside 
the cutoff radii, and the cutoff radii chosen in this scheme 
are normally larger than those in other pseudopotential 
generation schemes. The wave functions in the surface 
calculation were expanded in plane waves with a cutoff 
Ecnt = 55 Ry. To check if the cutoff energy is large 
enough for hydrogen, we calculated the total energy of a 
hydrogen molecule as a function of the cutoff energy. The 
difference in the total energy of II2 between i^cut = 55 Ry 
and 100 Ry was below 10 meV per H. For diamond the 
lattice constant and the bulk modulus were found to be 
3.53 A and 4.56 Mbar, respectively, in good agreement 
with the experimental values of 3.567 A and 4.43 Mbar. 
In the surface calculation, we used the theoretical lattice 
constant, which gave an ideal bond length of 1.53 A. 

Our slab model for the C(IOO) surfaces had the inver- 
sion symmetry and consisted of ten carbon layers and 
a vacuum region equal to six atomic layers (6.18 A). 
We varied the thickness of the crystal slab and the vac- 
uum region for the symmetric ( 1 x 1 ) : 211 phase , and found 
reasonable convergence of energy and geometry for our 
supercell. When varying the vacuum thickness, maxi- 
mum changes were less than 0.02 eV/(l x 1) in energy and 
about 0.5 % in bond length and angle using i^cut = 55 
Ry. The atomic positions at the two innermost layers 
were kept fixed. For the Brillouin-zone integration we 
used 6x6x2,3x6x2, and 2x6x2 grids in the Monkhorst- 
Pack special point scheme for the 1x1, 2x1, and 3x1 
phases, respectively. For the metallic (2xl):1.5II phase, 
a 4x8x2 grid was also tested to get better convergence. 
We determined the hydrogen and relaxed carbon posi- 
tions by minimizing the energy and Hellmann-Feynman 
forces. In our calculations, we used the force convergence 
criterion of 0.01 Ry/a.u., which gave a reasonable con- 
vergence in energy and geometry. Using a smaller force 
convergence criterion of 0.005 Ry/a.u. in the calculation 
of the symmetric (lxl):2II phase gave only a maximal 
change of 0.5% in the final bond length and an energy 
difference of 0.02 eV per (1x1) surface cell. 

For the bare surface and the (2x 1):1.5II phase, we used 
Gaussian broadening with a width of 0.1 eV to accelerate 
the convergence in the fc-point sum, because these struc- 
tures have dangling bonds on C atoms in intermediate 
or final configurations. We used a precondition Lanc- 
zos subspace-diagonalization procedure to solve for the 
eigenvalues and eigenvectors of the Hamiltonian matrix 
for a given LDA potential, and a modified quasi-Newton 
scheme to find the self-consistent potential. 
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III. RESULTS AND DISCUSSIONS 
A. Structures 

We have studied structures of different coverages and 
atomic configurations for the hydrogen-covered C(fOO) 
surfaces, including: (a) the bare reconstructed 2x1 sur- 
face, (b) the (2 X 1):H surface, (c) the (2xl):1.5H surface, 
(d) the (3xl):1.33H surface with alternating monohy- 
dride dimer and dihydride units, and (e) the symmetric, 
(f) canted, and (g) rotated (lxl):2H surfaces. We note 
that the unsaturated (2xl):1.5H surface was identical 
to the (lxl):1.5H surface considered by Furthmiiller et 
al}"^. We used the (2x1) notation to represent the peri- 
odicity of the surface unit cell. The top and side views 
of these structures are shown in Fig. 1. 

The calculated bond lengths for the top layers are 
given in Table I and compared with previous calcula- 
tions. Our C— H and top-layer C— C bond-length results 
were almost identical to those of earlier LDA calculations 
of Furthmiiller et al}"^ (using ultrasoft pseudopotentials 
and the LDA) and Zhang et al}'' (using ab tmtw MD) 
for structures they had considered. 

For the bare 2x1 surface, we obtained the symmet- 
ric dimer as the ground state. We started with a buck- 
led structure and the system relaxed to the symmetric 
one. This symmetric dimer structure is in contrast to 
that of Si(lOO) with a c(4x2) reconstruction composed 
of asymmetric dimers.^^ The origin of asymmetry has 
been identified as a result of competition between the 
TT bond and Jahn- Teller-like distortion. ^'^'^^ For the bare 
C(100)(2x 1), the filled 7r-bonding states were found to be 
separated from the empty antibonding tt* states even for 
symmetric dimers. Therefore, no Jahn- Teller type dis- 
tortion occurred. Similar to previous findings, the C— C 
dimer bond length on the bare C(100)(2x 1) surface (1.37 
A) was close to that of a typical carbon double bond. 
This strong tt and a bonding disappeared upon hydro- 
gen adsorption in the (2xl):H phase in which the C— H 
bond was formed. The C— C dimer bond length signifi- 
cantly lengthened to 1.63 A, being even larger than the 
bulk C— C bond length. 

For the (3xl):1.33H phase the bond lengths of the 
monohydride and dihydride were very similar to corre- 
sponding ones in the (2xl):H and symmetric (lxl):2H 
phases, respectively, as shown in Table L This structure 
is expected to have a low energy because the large steric 
repulsion caused by hydrogens in the (lxl):2H phase is 
greatly reduced by the presence of the dimers.^® The 
energetics will be considered in the next section. The 
(2xl):1.5H phase was considered by Furthmiiller et al}"^ 
as a possible high-coverage phase. We obtained almost 
identical geometries compared with their results. Since 
unsaturated dangling bonds exist in this structure, it is 
not energetically favorable as will discussed in the next 
section. 



Top view Side view 

[Oil] [100] 




FIG. 1. Top and side views of C(IOO) surfaces with dif- 
ferent hydrogen configurations: (a) bare 2x1, (b) (2xl):H, 
(c) (2xl):1.5H, (d) (3xl):1.33H, (e) symmetric (lxl):2H, (f) 
canted (lxl):2H, (g) rotated (lxl):2H. The bond lengths, in 
units of A, and bond angles are shown. Comparisons with 
previous theoretical results are given in Tables I and II. 

For completeness we also considered the structures of 
the (lxl):2H surfaces. A symmetric configuration [Fig. 
1 (e)] has a large steric repulsion between two hydrogens 
attached to adjacent carbon atoms. This repulsion can 
be reduced by a canted-row structure obtained from the 
symmetric dihydride by rotating the CH2 groups around 
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TABLE I. Calculated bond lengths (in A) in top two layers for various hydrogen-covered C(IOO) surfaces compared with 
results from other density-functional calculations. More calculated bond lengths are shown in Fig. 1. For comparisons with 
other previous theoretical studies, see Ref. 14. 





This work 


ReL 14"^ 




ReL 15 


ReL 16^ 


Lattice constant 


3.53 


3.53 


3.56 


3.62 


Bare 2x1 


symmetric 


symmetric 


symmetric 


buckled 


(iii(C-C) 


1.37 


1.37 


1.38 


1.40 


rfl2(C-C) 


1.50 


1.50 




1.52, 1.55 


(2xl):H 










(i(C-H) 


1.11 


1.10 


1.11 


1.17 


dii (C — C) dimer 


1.63 


1.61 


1.63 


1.67 


rfl2(C-C) 


1.53 


1.53 


— 


1.59 


(3xl):1.33H 










(i(C— H) monohydride 


1.10 


— 


— 


1.18 


(i(C-H) dihydride 


1.10 


— 


— 


1.17 


dii (C — C) dimer 


1.58 






1.62 


rfi2(C — C) monohydride 


1.53 






1.58 


7 / /' i x 1 \ 1 • 1 1*1 

rfi2(C — C) dihydride 


1.51 






1.60 


(2xl):1.5H 










(i(C— H) monohydride 


1.09 


1.09 


— 


— 


(i(C-H) dihydride 


1.11 


1.11 






rfi2(C — C) monohydride 


1.48 


1.47 






rfi2(C-C) dihydride 


1.50 


1.50 






(lxl):2H (symmetric) 










d(C-H) 


1.05 




1.06 




rfl2(C-C) 


1.50 








(lxl):2H (canted) 










d(C-H) 


1.11, 1.05 




1.04-1.13 




rfl2(C-C) 


1.51 









^Furthmiiller et al., self-consistent DFT 
•"Zhang et al., ah initio MD-DFT 
^Yang et al., non-self-consistent DFT 



an axis passing through the second layer carbon atoms 
[Fig. 1 (f)]. A rotation of ~10° increases the distance 
between hydrogen atoms in the same (neighboring) dihy- 
drides from L40 to L55 A (from LIO to L44 A), thereby 
lowering the energy. This canted dihydride structure was 
found to be 0.42 eV/(l x 1) lower in energy than the sym- 
metric dihydride. Alternatively the steric repulsion can 
be reduced by rotating the dihydride units around an axis 
perpendicular to the surface [Fig. 1 (g)], but we found 
that the energy gain was quite small, only 0.01 eV per 
1x1 unit cell. Also, we considered an alternating 1x2 
canted structure, which turned out to be 0.07 eV/(lxl) 
higher in energy than the original 1x1 canted one. 



Since different lattice constants were used in different 
theoretical calculations, a comparison of angles may be 
helpful. The angles H— C— C and H— C— H are given in 
units of degrees in Table II and shown in Fig. 1, along 
with the canting angle in the canted (lxl):2II phase. 
Our results are also similar to those of Furthmiiller et 
al.^^ and Zhang et al.^'' for the structures they had con- 
sidered. 



TABLE 11. Calculated angles in units of degrees (shown 
in Fig. 1) in various hydrogen-covered C(IOO) surfaces com- 
pared with results from some previous theoretical studies. 
The H— C— H angle for the canted dihydride phase of Ref. 
15 is for an alternating canted and twisted 1x2 structure. 





This work 


ReL 14 


ReL 15 


ReL 16 


(2xl):H 










H-C-C 


113 


113.3 


113.1 


113.9 


(3xl):1.33H 










H-C-H dihydride 


103 






108.7 


H-C-C monohydride 


110 






108.7 


(2xl):1.5H 










H-C-H 


105 


103.1 






(lxl):2H 










H-C-H (symmetric) 


83 




84.72 




H-C-H (canted) 


91 




93-96 




Canting angle 


11 
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B. Energetics 



Most previous calculations agreed that the dihydride 
structure had a high energy because of the steric repul- 
sion between two nearby hydrogen atoms, and that the 
(2xl):H monohydride structure, with all dangling bonds 
saturated, was a reasonably stable phase. Some calcu- 
lations also considered other competitive structures with 
higher hydrogen coverages. Furthmiiller et al}"^ stud- 
ied the (2xl):1.5H phase. Yang et al}^' using a nonself- 
consistent approach and Skokov et al}'' using a semiem- 
pirical Hamiltonian reported that the (3xl):1.33H struc- 
ture with alternating monohydride dimers and dihydride 
units was quite low in energy. Zheng and Smith^® found a 
strained symmetric dihydride configuration as stable us- 
ing slab-MINDO (modified intermediate neglect of differ- 
ential overlap). In order to compare stability among dif- 
ferent phases with different coverages, we will first adopt 
a chemical potential analysis^^ by considering the system 
in equilibrium with a hydrogen reservoir characterized by 
a chemical potential //jj- Later, the desorption energy re- 
sults will be discussed. 

The formation energy SI was defined as^^ 



Q = E + 



- ncnc - "H/"H, 



(1) 



where E was the total energy at zero temperature, Eq the 
zero-point energy, and nc and ripj the numbers of carbon 
and hydrogen atoms, respectively. Since the systems un- 
der study were assumed to be in chemical and thermal 
equilibrium with bulk carbon, we used the bulk energy 
as the chemical potential of carbon, jic- The zero-point 
energy included was from C— H vibrational modes. Since 
the hydrogen is always bonded to a single C atom, so it 
is reasonable to assume a zero-point energy proportional 
to the number of H atoms:^^ Eq = rifjeo. We estimated 
eo by the sum of vibrational frequencies of a typical hy- 
drocarbon CH4 and set it to 0.293 eV (0.0215 Ry) per 
hydrogen.^® Temperature dependence of the formation 
energy was mainly incorporated through the chemical po- 
tential of hydrogen. Although CVD is a highly nonequi- 
librium process, the study of equilibrium systems would 
be the first step in understanding stability of different 
phases during the growth. 

The calculated formation energies per 1x1 unit cell for 
different phase are shown in Fig. 2 as a function of the 
H chemical potential. The horizontal line corresponds 
to the bare 2x1 surface and the formation energy takes 
the value of the surface energy per 1x1 unit cell. Our 
result of 2.18 eV was close to the value of 2.12 eV ob- 
tained in Ref. 14. From Fig. 2, we can see that stability 
of each phase depends on the chemical potential of hydro- 
gen. As the H chemical potential increases, the bare 2x1, 
(2xl):H, (3xl):1.33H, and the canted dihydride phase 
successively become the most stable phase. The phase 
boundaries, indicated by arrows in Fig. 2, are located at 
— 17.19 eV, —15.66 eV, and —13.54 eV, respectively. This 
behavior is similar to that in the H-covered 
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FIG. 2. Formation energies Q [Eq. (1)], in units of eV per 
1x1 cell, for various phases as a function of hydrogen chem- 
ical potential jUh. As jUh increases, the bare 2x1, (2xl):H, 
(3xl):1.33H, and canted (1x1) phases successively become 
the phase with the the lowest formation energy. Arrows de- 
note the phase boundaries at —17.19 eV, —15.66 eV, and 
— 13.54 eV. The vertical dotted line at —15.49 eV, just be- 
side the cross point of the (2xl):H and (3xl):1.33H phases, 
indicates the chemical potential which gives a zero formation 
energy for CH4. See text for details. 

Si(lOO) surfaces reported in Ref. 27. The formation en- 
ergies for two selected chemical potential values, —0.957 
Ry (-13.02 eV) and -1.127 Ry (-15.33 eV), are listed 
in Table III, along with the total energies for each phase. 
The former is the energy of a hydrogen pseudoatom cal- 
culated with the local spin-density approximation (LSD), 
while the latter is the calculated energy per hydrogen in 
II2 with the zero-point energy of the vibrational mode 
(4400 cm~^)^^ included. In the first case, the canted 
(lxl):2H is most stable, while the (3xl):1.33H has the 
lowest formation energy in the second case. 

Not all the surface phases can be realized after taking 
into account other competing processes on the surface, 
such as the formation of hydrocarbons, within the same 
framework. We calculated the formation energy of CII4 
with the same theoretical methods using a supercell. The 
vertical dotted line in Fig. 2 near the crossing point of the 
(2xl):II and (3xl):1.33II phases represents a hydrogen 
chemical potential at which CII4 molecules can be formed 
with no cost in energy. The value is [£'(CIl4)-|-4eo— /Uc]/4 
= —15.49 eV. On the right hand side of the vertical line, 
the formation energy of CII4 is negative, namely, CII4 
forms spontaneously from the C(IOO) surfaces in the pres- 
ence of the hydrogen reservoir. In contrast, the left side 
of the vertical line has a positive formation energy for 
CII4. The most stable surface phases for this region are: 
the bare surface when //jj is smaller than —17.19 eV, the 
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TABLE III. Total energy E (in Ry) and formation energy Q (in eV, for two selected values of the hydrogen chemical 
potentials) per (1x1) unit cell for various C(IOO) surfaces. 





Bare 2x1 


(2xl):H 


(3xl):1.33H 


(2xl):1.5H 


(lxl):2H 
symmetric 


canted 


E 


-57.029 


-58.314 


-58.705 


-58.824 


-59.352 


-59.383 




2.18 


-2.00 


-2.88 


-2.28 


-2.81 


-3.23 




2.18 


0.31 


0.20 


1.19 


1.81 


1.40 



^using jUH = —13.02 eV, the LSD value for a free hydrogen atom 
•"using = |£(H2) = -15.33 eV 



TABLE IV. Desorption energy, in units of eV, of a hydrogen atom (see text) for various hydrogen-covered surfaces. The 
value for the (2xl):H phase is in reference to the bare (2x1) surface, while the values for other phases are with respect to the 
(2xl):H phase. The listed values were obtained using the LSD atomic energy of —13.02 eV, and did not include the zero-point 
energy of C— H vibrations. 





(2xl):H 


(3xl):1.33H 


(2xl):1.5H 


(lxl):2H 

symmetric 


canted 


This work 


4.47 


2.93 


0.86 


1.10 


1.52 


ReL 14 


4.54 




1.02 






ReL 15 


5.49 






1.13 




ReL 16 


6.20 


3.90 









(2xl):H phase when —17.19 < /Uh < —15.66 eV, and a 
small window for the (3xl):1.33H phase when —15.66 < 
//fj < —15.49 eV (see Fig. 2). The dihydride phase re- 
gion is totally excluded, as well as most of the 3x1 phase. 
Thus for a large //jj range (about 1.5 eV) the (2xl):H 
phase has the lowest formation energy, while the range for 
the (3xl):1.33H phase is quite limited (about 0.17 eV). 
In comparison, for the hydrogen covered Si(lOO) surfaces 
the hydrogen chemical potential that gives zero formation 
energy for SiH4 lies within a region where the dihydride 
is most stable. 

In general, the formation energy of defects involving 
higher hydrides will be reduced as //jj increases. There- 
fore, the (lxl):2II and (3xl):1.33II phases, which would 
become stable at larger //jj, are expected to be more 
susceptible to defect formation than the (2xl):II phase. 
These defects increase disorder on the surface, and may 
yield a 1 x 1 LEED pattern due to the underlying lattice 
as viewed through the disordered top layer. The rela- 
tively small //fj region for a stable 3x1 phase and the 
consideration of defect formation might explain why the 
C(100)(3xl):1.33H has not been observed so far.^s 

To answer the question whether other nxl phases can 
be realized, we estimated their energies by a simple model 
used in Ref. 27. It was based on the assumption that 
the formation energy could be approximated as a sum 
of the formation energies of isolated structural subunits 
and the interaction energy between these subunits. De- 
fine U as the difference in formation energy between the 
canted 1x1 phase and an isolated dihydride row unit. It 
is the energy arising from interactions between neighbor- 
ing rows of dihydrides, including the relaxation energy 
due to the CII2 rotation. We estimated that U = 1.41 



eV, compared with that of 0.10 eV for Si(100).27 This 
large repulsive interaction resulted in the wide range of 
the 3x1 phase. Using this model, the stability of any 
nxl structure consisting of a mixture of monohydride 
and dihydride units can be estimated. One type of them 
(n > 4) contains one monohydride row and n — 2 canted 
dihydride rows in each cell, and another consists of one 
dihydride row and (n — l)/2 monohydride rows. It can 
be easily shown that both of these nxl structures for 
C(IOO) surfaces are unstable relative to the phases we 
considered, same as concluded in Ref. 27 for Si(lOO). 

Another way to examine the energetics is to calculate 
the desorption energy for these systems. Following pre- 
vious theoretical studies, the desorption (adsorption) en- 
ergy AEn per hydrogen atom from the (2 x 1) :II surface is 
defined as the energy difference per surface site between 
the hydrogenated, reconstructed surface and the bare, re- 
constructed surface plus a free hydrogen atom We used 
the LSD energy of —13.02 eV for a free hydrogen atom. 
The desorption energies per hydrogen atom for higher- 
coverage structures were calculated with reference to the 
(2xl):II phase, instead of the bare surface. They can 
also be obtained by taking the difference of correspond- 
ing formation energies listed in Table III, and correct- 
ing for the zero-point energy and differences in coverage. 
Table IV summarizes these results and compares them 
with those obtained from other density-functional calcu- 
lations. More realistic desorption energy values should 
include the zero-point energy of C— H vibrational modes 
(about 0.29 eV per hydrogen) and reference to the energy 
of II2. Those values can be obtained by shifting down our 
results in Table IV by a constant of 2.60 eV. 

From Table IV, we see that our calculated desorption 
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TABLE V. Calculated vibrational energies of the H— C stretch mode in the (2xl):H phase. Results are obtained using the 
reduced-mass approximation and the masses of and ^^C. Energies are given in units of cm~^, while the values in parentheses 
are in meV. 





Theory 

This work 


Ref. 40 


Ref. 11 


Experiment 
Ref. 12 


Ref. 10 


Ground to first 


2835 (351) 


2750 


2915 


2910 


2928 


Ground to second 


5550 (688) 










Anharmonicity 


-120 (-15) 











energies are close to those of Furthiiller et al}"^ for sys- 
tems both considered. A larger difference was found 
between our (2xl):H result and that of Zhang et al.,^'' 
which may in part result from the F point sampling in 
Ref. 15. The difference between our results and those 
from the nonself-consistent calculation by Yang et al}^' 
was even bigger. 



C. Vibrational frequencies for the monohydride 
surface 

We have studied the hydrogen vibrational frequencies 
for the C(100)(2x 1):H surface within the frozen-phonon 
approximation by considering the energy as a function 
of hydrogen displacement from its equilibrium position 
along the stretch direction, with the carbon atoms fixed. 
This is the ordered structure observed experimentally 
on smooth surfaces so far. Since the hydrogen potential 




EIG. 3. Calculated energy for the H— C bond stretch mode 
in the C(100)(2 x 1):H phase as a function of H displacement r 
(in units of A) from its equilibrium position along the stretch 
direction. The solid line is a fourth-order polynomial fit, while 
the open circles represent calculated energies. The first three 
levels of vibrational energies are shown as horizontal lines. 



could not be well fitted by a harmonic form at larger 
atomic displacement, we fitted it by a quartic polyno- 
mial 

V(r) = Vo + LVlAr"^ - 2.310r^ + 2.292r'' (2) 

where V{r) is in Ry and the displacement r is in A. The 
plot is shown in Fig. 3. The Schrodinger equation of a 
hydrogen atom with this anharmonic potential was then 
solved numerically. The mass of the hydrogen atom must 
be reduced because of substate motion. In our calcula- 
tion, it was replaced by the reduced mass of a hydrogen 
atom and a carbon atom. It has been shown^^ that the 
effective-mass approximation is sufficient for the H— C 
stretch mode on the hydrogen covered C(lll) surface. 
Higher order corrections were estimated to be only 0.35 
% in Ref. 39 by a classical model with the spring constant 
between the top-layer carbon atoms and the substrate 
calculated from first principles. We performed a simi- 
lar calculation for the C(IOO) (2xl):II surface and found 
that the higher-order correction to be about 0.3 % within 
the simple classical model. The calculated energy sepa- 
rations between the ground and first excited states are 
shown in Table V, compared with other theoretical and 
experimental values. Our result was about 3 % smaller 
than the experimental values, similar to the difference 
between the LDA result and experiment for hydrogen 
on C(lll).^® The local-density approximation may be 
responsible for this difference, or the reduced-mass ap- 
proximation may have a larger error than estimated by 
the simple classical model. The calculated H— C stretch 
mode frequency is larger on C(IOO) than on C(lll) by 
about 3 %, same as the ratio of measured values for the 
two surfaces, although the calculated bond-length differ- 
ence on these two surfaces is no more than 1 %. We also 
examined the anharmonicity of the stretch mode, which 
is usually measured by the difference between the excita- 
tion energy from the ground state to the second excited 
state and twice the energy from the ground sate to the 
first excited state. The result, shown in Table V, was 
about 2 %, again similar to what was found for hydrogen 
on C(lll).39 

IV. SUMMARY 

We have studied the bare and hydrogen-covered C(IOO) 
surfaces using pseudopotentials within the LDA frame- 
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work. Different hydrogen coverages from one to two ML 
were considered. First, total energies and geometries of 
each phase were calculated. Relative stability between 
different phases with different coverages was investigated 
by examining the formation energy as a function of the 
hydrogen chemical potential. Taking into account the 
formation energy of CH4, two stable phases, (2xl):H 
and (3xl):1.33H, were found in different chemical po- 
tential ranges. The former was found in a wider range of 
— 17.19 < //H < —15.66 eV, while the latter in a relatively 
small region of —15.66 < //h < —15.49 eV. The region 
in which the (3xl):1.33H has the lowest formation en- 
ergy was originally much larger, but was significantly re- 
duced when the formation energy of CH4 was taken into 
account. This small region coupled with a higher pos- 
sibility of defect formation might explain why the 3x1 
phase has not been observed so far. Finally, we calcu- 
lated the vibrational energies of C— H stretch mode for 
the (2xl):H phase, and found them to be slightly larger 
than those on the C(lll) surface. 
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